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1.  INTRODUCTION 


The  small  wavelength  of  optical  radiation  (_  1/2  pm)  is  largely  responsible  for  the  high 
resolution  of  optical  imagery.  Even  a  small  lens  such  as  the  one  in  the  human  eye  images  very  well 
because  it  is  about  3000  wavelengths  in  diameter;  its  resolution,  according  to  the  laws  of  physics, 
is  1/3000  radian  or  about  0.02  deg.  Imaging  with  the  long  wavelength  used  in  radar  (roughly  l-3() 
cm)  is  much  more  difficult.  Consider  a  5-cm  imaging  radar,  for  which  the  wavelength  is  100,0{X) 
times  longer  than  its  optical  counterpart.  To  achieve  the  same  resolution  as  the  eye  the  antenna  size 
must  be  more  than  150  meters,  which  is  nearly  the  length  of  two  football  fields.  Even  more 
outlandish  is  the  5  km  antenna  required  to  achieve  the  resolving  power  of  a  small,  2"  optical 
telescope. 

Adaptive  beamforming  (ABF)  is  the  solution  to  the  unusual  design  that  achieves  these 
objectives.  With  ABF  an  air  traffic  controller’s  radar  could  utilize  target  echoes  from  several 

microwave  receivers  distributed  around  the  airport  and  deliver  high  resolution  images  -  not  merely 
the  usual  radar  blips  -  to  a  monitor  scope.  The  resolution  would  be  of  small-telescope  quality, 
good  enough  perhaps  to  see  if  the  wheels  were  down.  The  same  type  of  system  -  called  a  radio 

camera  -  aboard  an  airplane  would  provide  a  passing-scene  display,  in  the  cockpit,  of  the  terrain 
below  with  resolving  power  comparable  to  the  human  eye.  This  display  would  operate  in  all 
weather  and  night  and  ^y. 

Adaptive  signal  processing  is  the  means  by  which  ABF  enables  the  radio  camera  to  achieve 
this  extraordinary  resolution.  It  deduces  the  errors  in  the  locations  of  the  receivers  that  are 
dstributed  around  the  airport  or  on  the  air  frame  and  automatically  compensates  for  them  in  the 
image  processing.  The  information  is  obtained  from  the  echo-wavefront  measurements  made  by 
the  receivers.  In  short,  the  system  is  dynamically  self-calibrating. 

This  year's  work  concentrated  on  three  tasks.  The  first  was  to  develop  a  generalized  ABF 
theory  for  the  class  of  spatial  correlation  algorithms.  This  class  is  the  largest  and  appears  to  be  the 
most  generally  useful.  A  single,  general  algorithm  has  been  written  that  encompasses  all  known 
members  of  the  class.  When  the  tolerance  theory,  currently  being  developed,  is  finished,  a 
complete  engineering  design  theory  will  be  available  for  high  resolution  imaging  radar. 

The  second  task  extended  the  resolution  of  microwave  imaging  radar  from  a  1 -meter 
resolution  cell,  which  is  satisfactory  for  airplane  targets,  to  the  15  cm  needed  for  recognition  of 
ground  targets.  A  3-band,  dual  polarization  radar  was  built  with  1-GHz  bandwidth  for  this  task. 
Prominent  specular  echo  resonances  from  sides  and  floor  appear  in  images  of  small  pick-up  tmcks 

at  L-band  (10-cm  wavelength)  while  the  same  sides  and  floor  scatter  the  incident  energy  -  appear 

rough  -  at  K^-band  (1.6  cm).  This  difference  in  the  backscattering  characteristics  enhances  the 
target  recognition  potential  when  images  from  different  frequency  bands  are  combined. 

Enhanced  target  detection  sensitivity  and  target  recognition  was  the  third  thrust.  We 

improved  detection  sensitivity  by  the  order  of  20  dB  -  a  factor  of  100  -  with  an  electronic  3-D 
display  of  our  high  resolution  microwave  images  of  aircraft.  We  also  learned  that  simple  neural 
networks  can  recognize  the  images,  even  when  they  are  fragmented  and  noisy.  Our  current 
objective  is  to  automate  a  combination  of  the  two  techniques. 
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2 .  A  GENERALIZED  SPATIAL  CORRELATION-BASED  ALGORITHM 
FOR  SELF-COHERING  LARGE  ANTENNA  ARRAYS 


Many  potential  applications  exist  for  high-resolution  radar  such  as  direction 
finding,  high-accuracy  tracking,  target  counting,  and  high-resoludon  radar  imaging. 
All  of  these  applications  require  the  use  of  large,  thinned,  random  or  periodic 
antenna  arrays.  Many  uncertainties  exist  in  such  large  antenna  systems.  For 
example,  exact  element  positions  are  generally  not  known  because  of  surveying 
problems  or  flexing  of  the  large  antenna  structure.  In  addition,  different  channel 
lengths  and  impedance  mismatches  within  the  channels  will  cause  electrical 
distortion  in  the  array  system. 

Over  the  past  20  years  self-calibration  algorithms  have  been  developed  that 
require  nothing  more  than  data  acquired  from  a  target-of-opportunity  [1].  No  a 
priori  knowledge  concerning  the  environment  is  required.  As  this  work  developed, 
two  classes  of  algorithms  have  been  identified:  the  Dominant  Scatter  Class  and  the 
Spatial  Correlation  Class.  This  work  shows  that  the  entire  spatial  correlation  class 
of  algorithms  can  be  characterized  by  one  general  algorithm  called  the  Generalized 
Spatial  Correlation  Algorithm  (GSCA).  All  of  the  existing  correlation-based 
algorithms  are  shown  to  be  different  realizations  of  the  GSCA. 

In  this  paper  the  GSCA  is  derived.  Then  the  problem  of  calibrating  an  array 
with  this  algorithm  using  data  acquired  in  the  near  field  of  the  array  is  addressed. 
This  is  extremely  important  since  most  high-resolution  systems  operate  in  the  near 
field  of  the  array.  Finally,  the  relationships  between  the  existing  spatial  correlation 
algorithms  and  the  GSCA  tu'e  presented. 

Let  a  random  source  distribution  s(u)  be  illuminated  by  a  transmitting 
antenna  having  radiation  pattern  fT(u).  If  the  antenna  pattern  of  the  receiving  array 

is  fR(.u;  then  the  complex  image  of  the  source  distribution  is  §r(u)  =  fR(u)*fT(u)s(u) 

where  *  represents  convolution,  and  the  expected  image  intensity  is  E[Isr(u))I^].  If 
the  data  set  is  collected  over  K  rangebins  and  the  source  distribution  is  statistically 
homogeneous  over  these  rangebins  then  an  unbiased  estimate  of  the  expected  image 

K 

intensity  can  be  formed  as  g(u)  =  ~  1  IsR|^(u)I^  where  k  is  the  rangebin  index. 

k=l 

Define  p(u)  =  E[lfn(u)^f.j-(u)s(u)|2]  =  EflsoCu))!^]  to  be  the  expected  undistorted 

image  intensity.  .Now  consider  maximizing  the  function  J  =  |p(u)g(u)du.  The 

maximum  of  this  function  will  be  reached  when  g(u)=ap(u),  thus  driving  the  array 
towards  calibration.  However,  because  of  the  Fourier  relationship  between  an 
incoherent  source  distribution  in  the  far  field  and  the  spatial  correlation  of  the 
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received  radiation  in  the  array  aperture  J  =  F'lpgi^^Q  where  F'l  indicates  the 

inverse  Fourier  transform  and  L  is  the  discrete  lag  index  of  the  spatial  correlation 
function  measured  in  the  array  aperture.  It  can  be  shown  that  J  evaluates  to 

N-l 

J=ZwpU)^U)  (1) 

L^\ 

where  Wp(£)  is  a  weighting  function  determined  as  the  inverse  Fourier  transform  of 
p(u)  and  represents  an  esrimate  of  R(it)  (the  true  correlation  value)  obtained  by 
averaging  all  estimates  of  R(i)  available  in  the  array.  Equation  (1)  represents  the 
GSCA  objective  function,  maximization  of  which  drives  the  array  toward 
calibration.  Following  Hamaker  et.al.  [2]  maximization  of  J  enforces  a  stationarity 
constraint  on  the  correlation  values  of  lag  L  measured  at  all  available  pairs  of 
elements  in  the  array.  If  s(u)  is  a  wideband,  stationary  process  Wp(jd)  depends 
only  on  the  inverse  Fourier  transforms  of  the  transmit  and  receive  power  patterns. 


The  GSCA  objective  function  is  a  multiple  lag  objective  function  since  the 
sum  in  (1)  runs  between  L=\  and  jt=N-l.  The  GSCA  objective  function  has  been 
developed  under  the  assumption  of  far-field  conditions.  If  the  data  set  is  instead 
collected  from  the  near  field  then  the  data  .set  is  spatially  nonstationary  and  the 
theoretical  spatial  correlation  function  is  no  longer  independent  of  the  position  in  the 
array  at  which  the  correlation  is  calculated.  Consequently  there  is  no  reason  to 
enforce  a  stationarity  constraint  on  the  measured  correlation  values  of  various  lags. 
However,  depending  on  how  far  into  the  near  field  the  data  set  originates  certain 
low-order  correlation  lags  will  be  essentially  stationary  across  the  array  [3]. 
Therefore  the  GSCA  can  still  be  applied  to  calibrate  an  array  using  near-field  data. 
However,  potential  traps  do  exist.  The  weighting  function  Wp(;t)  is  based  solely  on 
far-field  assumptions  and  can  very  possibly  weight  heavily  certain  correlation  lags 
that  are  nonstationary  across  the  array  because  near-field  conditions  exist.  To 
account  for  this  potential  problem  a  .second  weighting  function  must  be  introduced 
into(l): 


N-l 

J  =  X  WNF(l)Wp(i)ft{;6) 


(2) 


L=\ 


The  weighting  function  Wf^jpU)  may  be  nothing  more  than  a  rectangular  window 

that  is  zero  for  lags  L  greater  than  the  highest-order  essentially  stationary  lag  and 
unity  for  all  lower-order  lags.  The  highe.st-order  es.sentially  stationary  lag  can  be 
determined  by  .solving  the  quadratic  inequality  l4j 


(3) 


u.- 


ov 


cJ^ 


for  L.  where  N  —  the  number  of  array  elements,  X  =  the  wavelength,  R  =  the 
distance  to  which  the  array  is  focu.sed,  d  =  the  interelement  distance,  and  Ufov  =  the 
extent  of  the  full  field  of  view  of  the  transmitter  beam.  The  solution  of  the  strict 
equality  of  (3)  yields  the  largest  lag  that  is  essentially  stationary  across  the  array.  It 
can  be  shown  that  the  phase  of  the  spatial  correlation  function  measured  in  the  array 
from  radiation  originating  in  the  near  field  can  be  .separated  into  two  terms;  one  that 
depends  only  on  the  lag  L,  and  another  (called  the  nonstationary  term)  that  depends 
on  both  the  lag  and  the  position  in  the  array  .  A  correlation  lag  is  considered  to  be 
essentially  stationary  across  the  array  when  the  nonstationary  term  is  less  than  a 
prescribed  value.  It  has  been  found  experimentally  [4]  that  a  nonstationary  term 
less  than  rt/io  radians  has  a  negligible  effect  on  the  calibration  process.  Using  n/io 
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radians  yields  the  (0. 15)  term  in  (3).  This  term  will  vary  depending  on  the  upper 
bound  chosen  for  the  nonstationtiry  phase  term  of  the  correlation  function. 

The  unit-lag  SCA  objective  function  developed  by  Attia  [5]  and  Subbaram's 
multiple-lag  SCA  objective  function  [3j  are  both  special  cases  of  the  GSCA's 
general  objective  function  (2).  The  objective  function  of  the  multiple-lag  SCA  can 

be  obtained  by  setting  wpU)  =  R(ji)  and  then  estimating  R(jt)  as  6.(Z)  since  R(i)  is 
unknown  [6].  The  unit-lag  SCA  objective  function  is  obtained  in  the  same  manner 
as  the  multiple-lag  SCA  but  the  near-field  weighting  function  wnfU)  is  zero  for  all 
lags  but  the  unit  lag. 

Insight  may  now  be  gained  by  considering  these  spatial  correlation-based 
algorithms  as  special  cases  of  the  GSCA.  For  example,  from  the  preceding 
discussion  on  using  the  GSCA  in  the  near  field,  it  should  be  expected  that  any 
spatial  correlation-based  algorithm  will  breakdown  when  calibrating  too  far  into  the 
near  field.  However,  as  the  calibration  moves  further  into  the  near  field,  wnf(1) 
will  be  unity  for  fewer  lags  until  at  some  point  it  will  be  unity  for  only  the  unit-lag. 
Moving  any  further  into  the  near  field  will  cause  all  lags  to  be  nonstationary  and  the 
GSCA  will  fail  completely.  Consequently  the  unit-lag  SCA  is  the  most  robust  (in 
terms  of  near-field  calibration)  spatial  correlation-based  algorithm.  Rarely  is 
calibration  done  so  tar  into  the  near  field  that  the  unit-lag  becomes  nonstationary, 
and  hence  the  success  Attia  achieved  with  it. 

The  GSCA  can  be  used  to  gain  even  more  dramatic  insight  into  Tsao's 
Energy  Conservation  Algorithm  (ECA)  (7l  which  has  previously  been  thought  to 
represent  another  class  of  self-calibration  algorithms.  The  algorithm  is  based  on 
Parseval's  Theorem  and  the  objective  function  used  by  Tsao  minimizes  the 
difference  between  the  total  energy  in  the  array  aperture  and  the  total  energy  in  the 
image  within  a  specific  angular  sector  determined  by  the  field  of  view  of  the 
transmitter  beam.  Since  the  total  energy  within  the  aperture  is  constant,  minimizing 
Tsao's  objective  function  is  equivalent  to  maximizing  the  total  energy  within  the 
specific  angular  sector  in  the  image.  Therefore  an  equivalent  ECA  objective 
function  is  [  1 1 


^  k=i 


'Au/2 


-iiu/2 


(u)F  du 


(4) 


However,  this  objective  function  is  nothing  more  than  J  with  p(u)  =  rect(Au). 
Therefore  the  ECA  i^  another  spatial  coirelation-based  algorithm  with  Wp(jt)  equal 
to  a  sine  function. 

Figures  1  and  2  demonstrate  the  effect  of  the  weighting  function  Wp(Z)  in 
(1).  The  array  was  calibrated  on  a  simulated  source  distribution  which  contained 
10  rangebins  each  with  2  point  sources  randomly  distributed  over  an  angular  extent 
of  0.167  radians  except  the  first  rangebin  which  had  one  source  located  at  zero 
radians.  Figure  1  shows  the  image  of  the  first  rangebin  obtained  after  calibrating 
the  array  with  the  GSCA  using  the  appropriate  weighting  function  wp(Z) 
determined  by  p(u).  This  diffraction  limited  image  clearly  demonstrates  the  array 
has  been  calibrated  correctly.  Figure  2  shows  the  image  of  the  same  point  source 
obtained  by  calibrating  with  the  GSCA  using  an  incorrect  weighting  function 
corresponding  to  a  rectangular  function  in  the  u  domain  of  much  less  angular  extent 
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than  is  illuminated  by  the  transmitter.  The  effect  of  using  the  incorrect  weighting 
function  is  seen  by  the  distorted  image  of  the  point  source  in  Figure  2  .  The  ECA 
was  developed  by  Tsao  as  an  independent  algorithm.  Now  with  the  understanding 
that  the  ECA  is  really  a  correlation  algorithm  it  can  be  seen  that  the  ECA  will  only 
work  well  in  simulations.  This  is  because  the  rectangular  function  p(u)  implies  that 
no  energy  was  radiated  or  received  outside  of  the  mainlobe  of  the  system.  This  is 
not  true  in  practical  systems  where  energy  is  always  radiated  and  received  outside 
of  the  mainlobe  by  sidelobes.  Consequently,  wp(jt)  equal  to  a  sine  function 
describes  an  unrealizable  situation. 


Figure  1.  Image  of  a  point  source  Figure  2.  Image  of  a  point  source 
located  at  u=()  with  p(u)  =  located  at  u=0  with  p(u)  = 

rect(().167  radians).  rectlO.0 1724  radians). 

The  GSCA  is  an  important  generalization  of  the  existing  spatial  correlation- 
based  algorithms.  The  generalized  objective  function  provides  a  basis  from  which 
analysis  applicable  to  the  entire  spatial  correlation  class  can  be  carried  out.  One 
such  important  analysis  is  the  derivation  of  a  tolerance  theory  for  this  class  of 
algorithms  similar  to  the  theory  that  already  exists  for  the  Dominant  Scatterer  class 
[8]. 
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3.  2-D  HIGH  RESOLUTION  MICROWAVE  IMAGING  OF 
GROUND  VEHICLES  WITH  A  UWB  RADAR  SYSTEM 


INTRODUCTION 


During  1989  a  system  was  developed  at  VFRC  for  the  detection  of  inert  military  landmines.^ 
Since  such  mines  are  very  Umited  in  size,  it  was  important  to  develop  resolution  cells  that  were  as  small  as 
possible.  An  S-band  radar  set  with  1  GHz  bandwidth  was  developed  for  this  purpose.  This  UWB  radar 
system  provided  a  range  resolution  of  15  cm  in  air.  It  was  only  a  short  time  after  the  completion  of  the 
mine  detection  program  when  this  equipment  was  applied  to  the  imaging  of  larger,  ground  based  surface 
targets. 


This  paper  will  explain  both  the  design  of  this  system  and  the  extensions  which  were  made  during 
1991  to  allow  operation  in  three  diverse  frequency  bands.  A  variety  of  target  images  from  the  resulting 
experiment  program  will  also  be  presented. 


IMPORTANT  FEATURES  OF  THIS  PROGRAM 

The  original  imaging  system  developed  in  the  mine  detection  effort  used  a  1  GHz  bandwidth  with  a 
center  frequency  of  3.5  GHz.  The  fractional  bandwidth  (0.286)  therefore  classifies  this  system  as  a  stand¬ 
alone,  ultra- wideband  radar. 

It  was  very  clear,  however,  that  the  wideband  video  and  digital  components  operated 
independently  of  the  center  microwave  frequency.  Equipment  from  a  variety  of  past  programs  was 
available  which  allowed  additional  front  ends  to  be  constructed  at  center  frequencies  of  9.6  and  16.8  GHz. 
Thus  the  same  targets  could  be  imaged  at  three  different  frequencies  representing  a  total  span  of  more  than 
two  octaves.  Although  1  GHz  bandwidth  systems  at  the  higher  two  frequencies  are  not  classified 
inherently  as  UWB  radars,  their  combination  forms  a  system  of  enormous  bandwidth  (although  measur^ 
at  discrete  intervals).  This  represents  a  concept  of  spectrum  conservation  which  may  prove  important  in 
the  future  growth  of  UWB  system  concepts. 

A  further  important  result  of  this  work  was  the  production  of  images  with  two-dimensions  as 
contrasted  with  the  traditional  range/amplitude  profiles.  As  will  be  shown  later,  one-dimensional  range 
scans  are  useful  and  provide  important  input  to  target  identification  and  classification  algorithms.  The  two- 
dimensional  image,  however,  adds  angular  resolution  to  provide  another  important  source  of  information. 


SYSTEM  DESIGN 

The  simplified  block  diagram  of  the  three-band  UWB  radar  system  is  presented  in  figure  3.  This 
design  insured  that  the  video  and  digital  sections  were  frequency  independent  in  order  to  minimize  the  cost 
of  system  development.  Only  one  modular  radar  system  was  constructed.  The  array  was  formed  by 
physically  moving  the  front  end  to  each  specified  sampling  position.  Each  of  the  three  systems  was  fully 
coherent,  and  the  experiments  were  virtually  time  insensitive.  This  design  also  permitted  great  flexibility 
in  reconfiguring  the  equipment  for  a  variety  of  experiments.  The  only  disadvantage  of  the  synthetically 
formed  array  was  that  only  stationary  targets  could  be  imaged.  This  did  not  prove  to  be  a  liability, 
however,  since  a  large  variety  of  suitable  objects  was  available. 

Although  three  frequencies  were  in  use,  the  block  diagram  for  each  band  was  the  same.  A  single 
microwave  oscillator  was  used  to  provide  both  the  transmitted  pulse  and  the  coherent  reference  for  the 
phase  detector.  The  1  nsec  pulse  (3  dB)  was  formed  by  a  New  England  Microwave  series  6500  pin  diode 
switch.  The  peak  power  level  from  this  device  was  on  the  order  of  1  mW.  A  free-running  HP  8082  pulse 
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Figure  3.  Block  Diagram  of  the  Three  Frequency  UWB  Radar  System. 


(a)  (b) 

igure  4.  Pickup  Truck  Target  (length,  5.4m,  width  2m,  length  of  bed  2.5m) 

(a)  Overhead  plan  view 

(b)  In  perspective  as  viewed  from  the  array 


generator  was  used  to  provide  the  gate  to  the  switch.  It  also  established  the  PRF  of  the  system  at  a  rate  of 
approximately  100  kHz. 

Amplification  of  the  power  level  was  accomplished  by  external  solid  state  or  TWT  amplifiers. 
Maximum  power  output  was  limited  to  10  Watts. 

Separate  hom  antennas  were  used  for  the  transmitter  and  receiver  in  each  systetiL  The  use  of  a  T/R 
switch  was  not  considered  due  to  the  very  close  range  of  the  targets.  The  gains  of  the  horns  for  the  three 
band  system  varied  from  18  to  24  dB.  Due  to  the  extreme  near  field  conditions  of  the  experiments,  the 
antennas  were  continually  repointed  at  the  target  complex  as  they  were  moved  to  each  successive  sampling 
location.  The  noise  figure  of  the  system  was  established  by  a  low  noise  amplifier  connected  to  the  output 
of  the  receive  hom. 

The  phase  detector  is  based  on  a  pair  of  dc -coupled,  wideband  mixers.  Two  circuits  were  built. 
One  for  S-band  alone  and  the  other  for  both  X  and  Ku-bands.  The  90  degree  line  stretcher  was  adjusted 
each  time  the  system  was  changed  between  the  two  higher  frequencies.  Since  the  inputs  to  the  phase 
comparator  were  of  identical  frequency,  the  baseband  output  was  a  dc  level  proportional  to  the  amplitude 
and  phase  of  the  input  signals.  This  level  persisted  for  the  length  of  the  radar  echo.  (1  nsec)  Thus  the 
filtered  output  of  the  phase  comparator  contained  the  quadrature  signals  required  by  the  image  processing 
software. 


The  transmitter  is  obviously  a  key  component  in  the  development  of  a  UWB  radar  system,  but  an 
effort  of  similar  magnitude  exists  in  providing  an  A/D  converter  capable  of  digitizing  waveforms  with  very 
wide  bandwidths.  The  HP  54120T  is  a  sampling  oscilloscope  that  is  operational  from  dc  to  20  GHz. 
Thus  not  only  can  the  actual  transmitted  RP  waveform  be  displayed,  but  the  baseband  I  and  Q  signals  can 
be  digitized  and  stored  in  RAM  for  downloading  to  the  control  computer.  The  sampling  interval  can  be  set 
as  low  as  10  psec,  thus  removing  the  speed  of  the  digitizer  as  a  limiting  factor  in  an  experimental  program. 
Waveforms  can  also  be  averaged  to  improve  the  SNR.  The  only  disadvantage  of  the  instrument  is  the  time 
required  to  complete  the  sampling  procedure  (approximately  10  seconds  for  x32  averaging).  Since  the 
targets  are  stationary,  however,  this  restriction  is  not  significant.  The  oscilloscope  is  connected  to  an 
HP98 16  desktop  computer  through  the  HPEB  interface  line.  This  allows  convenient  transfer  of  data  sets  to 
floppy  discs  for  permanent  storage. 

As  a  summary  of  the  operating  features  of  the  system  the  parameters  of  the  S-band  set  are 
summarized  below: 


center  frequency 
wavelength 
pulse  width  (3  dB) 
system  bandwidth  (3  dB) 
peak  transmitted  power 
noise  figure 
antenna  gain 
antenna  beamwidth 
minimum  range 
maximum  range 
array  type 


3.5  GHz 
0.0857  meters 
1  nsec 
1  GHz 

0.1  W  to  10  W  as  needed 

3dB 

18  dB 

20  degrees 

no  limit 

approximately  250  meters 
monostatic  synthetic  aperture 


EXPERIMENT  PROGRAM 

The  overall  objective  of  this  phase  of  the  program  was  to  study  the  images  of  complex  targets 
under  ideal  conditions.  In  previous  experiment  programs  the  targets  under  study  were  often  placed  in 
complex  surroundings  at  long  range.  In  order  to  maximize  the  image  quality  the  current  tests  were 
conducted  at  close  range  on  a  clear  section  of  field.  After  these  studies  have  been  completed,  it  will  be 
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very  important  to  repeat  them  in  more  typical  environments  with  both  foliage  and  ground  clutter.  The 
targets  were  positioned  approximately  10  to  15  meters  from  the  array,  and  the  angular  resolution  required 
to  match  the  15  cm  range  cell  was  on  the  order  of  10  mrad.  Each  anay  was  50  wavelengths  long,  and 
element  sampling  was  conducted  each  half  wavelength.  Thus  the  array  lengths  were:  4.24  m  at  S-band, 
1.55  m  at  X-band,  and  0.88  m  at  Ku-band. 

Since  the  system  had  already  been  constrained  to  image  only  stationary  targets,  the  synthetic 
apenure  technique  was  used  to  measure  sufficient  element  positions  to  form  a  filled,  periodic  array. 
Future  experiments  with  array  thinning  may  show  the  usefulness  of  this  time  saving  procodure. 

Both  transmitter  and  receiver  were  moved  to  each  successive  sampling  position  as  the  synthetic 
apermre  was  formed.  The  resulting  monostatic  array,  when  coupled  with  15  cm  range  resolution, 
provided  the  potential  for  very  detailed  images. 


RANGE/AMPLITUDE  EXPERIMENTAL  RESULTS 

Range/amplitude  presentations  have  traditionally  been  used  to  describe  radar  cross  section  profiles 
at  a  specified  aspect  angle.  This  is  a  very  useful  method  for  gaining  insight  into  certain  characteristics  of 
the  targets  under  study.  It  has  the  added  advantage  that  no  signal  processing  is  required,  thus  allowing  the 
results  to  be  observed  immediately.  In  fact,  research  at  VFRC  on  using  neural  networks  to  classify  targets 
from  range/amplitude  profiles  has  shown  promise  of  successful  operation. 

A  typical  target  used  in  this  program  was  a  conventional,  full  size  pick-up  truck.  The  following 
discussion  will  focus  on  this  object  in  order  to  provide  a  detailed  analysis  of  the  resulting  images.  The 
truck  was  initially  positioned  in  a  clear  section  of  field  on  the  VFRC  test  range  facing  away  from  the  array 
as  shown  in  figure  4. 

The  radar  equipment  viewed  the  target  area  from  a  platform  approximately  12  degrees  higher  in 
elevation.  S-band  was  the  first  frequency  chosen  for  the  tests.  The  equipment  was  configured  for  vertical 
polarization  and  positioned  at  the  array  center.  The  range  profile  shown  in  figure  5  was  observed.  Notice 
that  six  distinct  scattering  centers  are  observed  in  the  region  between  6  and  20  meters  in  range.  The  first, 
located  at  approximately  7  meters,  is  a  small  ( 1  foot  in  dimension)  comer  reflector  used  as  a  reference 
target.  Large  specular  reflections  from  both  the  tailgate  and  the  rear  of  the  cab  are  found  at  9  and  12 
meters,  respectively.  The  responses  between  14  and  20  meters,  however,  did  not  correspond  to  any 
targets  in  position  on  the  field.  It  was  soon  determined  that  they  were  from  multiple  reflections  between 
the  tailgate  and  the  cab  in  the  bed  of  the  truck. 

In  order  to  verify  that  this  was  the  correct  interpretation  of  these  echoes,  the  tailgate  was  lowered 
without  moving  either  the  truck  or  the  radar  set  This  eliminated  the  multiple  reflections  as  shown  in  figure 
6.  The  response  of  the  rear  of  the  cab  was  much  larger  in  this  case  since  the  scattering  surface  of  the 
tailgate  had  been  significantly  decreased. 

-After  replacing  the  tailgate  in  its  original  closed  position,  the  first  synthetic  apermre  experiment  was 
conducted.  During  this  procedure,  the  radar  set  was  moved  over  a  distance  of  4.24  meters.  This 
corresponds  to  a  variation  in  aspect  angle  of  approximately  25  degrees  at  the  rear  of  the  truck.  The 
multiple  reflections  were  observed  (with  varying  amplitude)  over  this  entire  angular  extent.  Variation  in 
elevation  was  tested  between  7  and  30  degrees.  Again,  the  multiple  reflections  were  identified  in  each 
experiment.  Their  amplimdes  peaked,  however,  between  10  and  15  degrees. 

Experiments  were  also  conducted  with  horizontal  polarization.  The  character  of  the  truck  response 
was  very  similar  to  that  observed  with  venical  polarization.  The  multiple  reflections,  however,  were 
significandy  decreased  in  level. 
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7.  X-band  rear  view. 


Figure  8.  Ku-band  rear  view. 


The  results  from  tests  conducted  at  X  and  Ku-bands  are  also  of  interesL  They  both  showed  very 
large  specular  reflections  from  the  tailgate  and  reduced  responses  from  the  cab.  At  these  higher 
frequencies,  the  multiple  reflections  were  so  small  as  to  be  almost  unidentifiable.  These  range/amplitude 
profiles  are  shown  in  figures  7  and  8  .  They  again  represent  the  rear  view  of  the  truck  from  the  array 
center. 


For  a  more  complete  understanding  of  the  scattering  properties  of  the  truck  it  was  examined  at  a 
variety  of  aspect  angles.  Shown  here  are  the  results  of  S-band  tests  where  the  truck  was  viewed  from  the 
rear  but  38  degrees  off  axis  (figure  9),  parallel  to  the  array  (figurelO,  and  from  the  front  (figurel  1).  In 
each  of  these  cases  at  least  one  multiple  reflection  was  clearly  identified.  It  should  be  noted,  however,  that 
large  differences  in  both  amplimde  and  multiple  reflection  delay  are  very  apparent  and  tend  to  verify  that 
the  one-dimensiond  radar  si^ature  of  the  vehicle  is  quite  dependent  on  aspect  angle.  Despite  their 
complexity,  these  high  resolution  amplitude/range  plots  are  nevertheless  very  irrformative  in  attempting  to 
classify  this  object  as  a  small  truck. 


TWO-DIMENSIONAL  IMAGING  PROGRAM 

Further  ability  to  identify  the  target  may  be  obtained  by  adding  a  second  dimension  to  the 
range/amplitude  display.  After  the  quadrature  target  echoes  have  been  measured  at  each  of  the  sampling 
positions  in  the  synthetic  aperture  array,  several  procedures  are  used  to  produce  a  two-dimensional  image. 
First,  the  ech^s  must  be  aligned  to  compensate  for  the  variation  in  range  as  the  radar  set  is  moved  along 
the  array.  This  is  strictly  accomplished  according  to  the  geometry  of  the  experiment.  Since  the  range  cell 
is  only  15  cm  in  depth,  however,  this  procedure  must  be  conducted  not  only  for  each  range  bin  but  also 
for  individual  scan  angles.  It  is,  therefore,  a  time  consuming  step  in  the  image  processing. 

The  phased  array  itself  is  produced  offline  using  adaptive  beamforming  to  compensate  for  errors  in 
the  position  of  the  radar  module  during  the  experiment.  This  is  a  blind-equalization  means  of  self¬ 
calibrating  a  distorted  antenna  array.  The  target  echoes  received  by  the  antenna  elements  have  been 
randomly  modulated  in  phase  by  the  element  position  errors.  The  adaptive  beamforming  processor 
extracts  information  from  the  measured  phase  of  the  reference  reflector  from  which  it  calculates  the 
required  correction  for  each  element's  measured  phase  values.  Figure  12  shows  the  point  source 
response  of  the  X-band  array.  Note  that  the  abscissa  in  this  plot  is  cross  range,  not  range  as  in  the 
previous  examples.  This  is  one  range  slice  of  the  80  which  are  used  in  forming  a  two-dimensional  image. 
It  covers  the  range  of  the  small  reference  comer  reflector,  and  demonstrates  the  10  mrad  angular  resolution 
of  the  array.  Attempts  to  process  this  data  set  without  using  adaptive  beamforming  resulted  in  a  distorted 
output  image.  This  verified  that  the  tolerance  on  element  location  was  not  adequate  during  the  experiment 
to  allow  for  the  use  of  conventional  array  processing  techniques. 

The  two-dimensional  images  for  the  pick-up  truck  viewed  from  the  rear  are  shown  in  figures  13 
(S-band),  14  (X-band),  and  15(Ku-band).  As  was  noted  in  the  range/amplitude  profiles,  the  multiple 
reflections  are  prominent  in  the  S-band  experiment  but  very  small  at  the  higher  frequencies.  The  reflection 
from  the  rear  of  the  cab  is  similarly  dependent  on  the  band.  As  the  frequency  rises,  its  amplitude 
decreases.  The  image,  in  fact,  becom  es  more  specular  and  less  recognizable  optically  as  the  frequency  in 
increased.  The  angular  dimension  in  these  plots,  however,  makes  visual  identification  of  the  target  as  a 
small  truck  possible. 

Another  advantage  of  the  two-dimensional  presentations  is  that  diversity  techniques  may  by  used  to 
further  improve  the  image.  Shown  in  figure  16  is  a  demonstration  of  combining  the  three  maps  shown 
previously  (figures  13,14,  and  1 5 )  to  form  a  composite  image  representing  the  characteristics  of  all  three 
frequencies.  The  multiple  reflections,  large  speculars,  and  target  detail  of  each  individual  frequency  range 
are  apparent  in  this  image.  In  effect,  this  plot  represents  a  UWB  radar  system  with  discretely  sampled 
bandwidth  of  over  14  GHz. 
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Figure  12.  X-band  specular  reflector 
image  illustrating  angular 
resolution  of  array  (10  mrad) 


Figure  13.  S-band  rear  view  image  in  two  dimensions.  (Pickup  truckl 
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Figure  16.  Rear  view  of  truck,  composite  image  created  by  summing  figures 
11,  12,  and  13.  The  i.mage  contains  the  important  features  of 
all  three  frequencies. 
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4.  ENHANCED  TARGET  DETECTION  USING 
STEREOSCOPIC  DISPLAY  TECHNIQUES 


INTRODUCnON 

In  an  experiment  performed  at  the  Valley  Forge  Research  Center,  it  was  observed  that 
target  detection  sensitivity  could  be  increased  by  using  high  resolution  microwave  imaging  in  two 
and  three  dimensions  fi  i].  Continued  research  during  this  period  comprised  three  tasks.  The  first 
was  a  critical  analysis  of  the  experimental  procedure  and  an  interpretation  of  the  results.  Secondly, 
a  smdy  was  made  of  three-dimensional  perception  and  display  in  order  to  gain  better  understanding 
of  the  mechanisms  involved.  Finally,  applications  were  considered  to  determine  where  this 
technology  may  prove  most  useful 

EXPERIMENTAL  ANALYSIS 

In  the  first  part  of  the  experiment,  detection  sensitivity  was  plotted  against  target-to-clutter 
ratios  (TCR)  for  0-dimensional  and  2-dimensional  targets.  The  definition  of  TCR  used  in  the 
experiments  is  given  as 


TCR  -  Power  per  Pixel  in  the  Target 

Average  Power  per  Pixel  in  the  Clutter 

The  reported  result  was  that  the  2-D  targets  could  be  detected  with  equal  probability  at  TCR's  1  IdB 
lower  than  the  0-D  targets.  In  the  second  part  of  the  experiment,  detection  sensitivity  was 
compared  for  3-D  and  2-D  targets.  This  resulted  in  equal  detection  at  TCR’s  an  additional  7dB 
lower. 

Correction 

During  this  period,  a  minor  error  that  bears  mentioning  was  discovered  in  the  first  part  of 
the  experiment.  In  the  0-D  case,  the  target  pixel  replaced  a  clutter  pixel  in  the  final  image. 
However,  it  should  have  been  added  coherently  to  the  clutter  as  in  both  the  2-D  and  3-D  cases.  As 
a  result,  the  target  power  used  to  determine  the  TCR  in  the  0-D  case  was  larger  than  it  should  have 
been.  We  can  determine  approximately  what  the  target  power  was  by  subtracting  from  this  value 
the  average  clutter  power.  The  correction  for  the  TCR  was  performed  and  the  new  0-D  curve  is 
displayed  in  Figure  17.  For  high  TCR's  the  correction  is  negligible  and,  therefore,  the  result 


reported  earlier  still  stands.  That  is,  for  high  TCR,  it  does  not  matter  significantly  whether  the 
target  replaces  or  is  added  to  the  clutter. 

Conclusions' 

In  the  earlier  report,  the  improvement  in  detection  sensitivity  was  attributed  to  the  increased 
target  dimensionality.  During  this  period,  we  determined  that  the  improvements  that  were 
wimessed  were  the  result  of  two  factors. 

In  the  first  part,  of  the  experiment,  the  improvement  is  not  merely  a  result  of  increased 
dimensionality,  but  of  the  increased  resolution  at  which  the  target  is  displayed.  This  is  easily  seen 
when  the  0-D  target  is  merely  caisidered  to  be  a  very  small  2-D  target.  Now  both  targets  have  the 
samp,  dimensionality,  but  the  target  with  the  finer  resolution,  the  2-D  target,  has  a  much  lower 
threshold  of  detection.  When  the  resolution  of  the  target  is  increased,  the  target  occupies  more 
pixels  in  the  image.  Since  each  pixel  conveys  additional  information,  it  becomes  easier  to  detect 
the  target  as  the  resolution  becomes  finer. 

The  0-D  case  can  then  be  interpreted  as  the  limiting  case  for  2-D  as  the  resolution  cell 
grows  larger  than  the  size  of  the  target  Therefore,  we  can  expect  that  the  2-D  curve  in  Figure  17 
will  approach  the  0-D  curve  as  the  resolution  cell  grows.  Similarly,  it  may  be  possible  to  shift  the 
2-D  curve  even  farther  away  from  the  0-D  curve  by  improving  the  resolution  of  our  imaging 
system. 

In  the  second  part  of  the  experiment,  the  improvement  is  a  result  of  the  observable  depth 
difference  between  the  target  and  the  clutter  in  the  three-dimensional  image  [12].  The  manner  in 
which  this  depth  difference  is  achieved  is  discussed  in  the  next  section.  For  now,  however,  we 
can  state  that  the  2-D  case  is  the  limiting  case  of  3-D  where  the  depth  difference  is  decreasing. 
Therefore,  it  is  important  to  determine  the  factors  necessary  to  produce  these  large  differences. 

TEIREE-DIMENSIONAL  PERCEPTION  AND  DISPLAY 

Because  of  the  obvious  benefit  of  three-dimensional  detection,  a  study  was  undertaken  to 
gain  a  better  understanding  of  the  mechanisms  involved  in  the  display  and  perception  of  three- 
dimensional  images.  What  follows  is  a  brief  introduction  to  some  key  points. 

Disparity 

When  we  look  with  two  eyes  at  an  object,  each  eye  sees  a  slightly  different  image  of  the 
object  since  the  viewing  positions  of  the  left  and  right  eye  are  different.  This  phenomenon  is 
known  as  stereo  vision  or  stereopsis  and  can  be  experienced,  for  example,  by  looking  at  your 
index  finger  while  alternately  opening  and  closing  each  eye.  What  you  will  see  is  that  the 
horizontal  position  of  your  finger  will  appear  to  shift  with  respect  to  the  background.  The  brain 
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uses  this  shift  to  judge  the  distance  between  your  finger  and  the  background.  The  mechanism  of 
judging  distance  based  on  the  perceived  differences  in  the  images  provided  by  each  eye  is  known 
as  binocular  parallax  or  binocular  disparity  and  is  further  demonstrated  using  Figure  18, 

When  the  two  eyes  focus  on  a  point,  denoted  by  A  in  this  figure,  the  eyes  rotate  until  the 
focused  point  is  projected  onto  the  centers  of  the  retinae,  labeled  A'.  All  other  points,  such  as 
point  B,  "fall  where  they  may"  in  the  retinal  image.  Notice  that  the  separation  of  A'  and  B'  is 
larger  in  the  right  retinal  image  than  in  the  left.  These  separations  are  typically  measured  as  angles 
and  are  denoted  by  ©i  and  ©2  in  the  figure.  The  disparity  of  a  point  is  defined  as  the  difference  of 
these  angles  and  is  give  notationally  as 


ti  =  ©i-©2 


It  can  be  shown  that  the  perceived  distance  between  the  points  A  and  B  is  directly  proportional  to 
the  disparity  of  point  B.  It  follows  that,  for  our  purposes,  larger  disparities  will  normally  result 
in  better  separation  between  the  target  and  background  clutter. 

The  disparities  that  the  brain  can  detect  vary  from  person  to  person  and  have  also  been 
found  to  be  dependent  on  the  size  of  the  stimulus  used.  Current  literature  suggests  that,  for  the 
large  targets  used  in  our  experiment,  depth  differences  should  typically  be  perceived  for  disparities 
ranging  from  about  10"  to  2°  [i3j.  For  disparities  above  this  range,  the  brain  experiences  difficulty 
fusing  the  images  to  form  a  three-dimensional  scene.  In  fact,  even  within  this  range,  fusion  may 
often  take  a  considerable  amount  of  time  for  large  disparities.  The  perceivable  range  of  disparity  is 
an  important  factor  to  consider  in  the  production  of  the  three-dimensional  images  displayed  in  our 
experiments. 

In  Figure  19,  we  demonstrate  how  the  3-dimensional  scene  of  Figure  1 8  can  be  reproduced 
by  showing  the  eyes  two  separate  2-dimensional  views  of  the  scene.  Notice  that  the  separations  of 
A’  and  B'  are  identical  to  those  experienced  in  Figure  18  .The  result  is  that  the  brain  is  "tricked" 
into  beheving  that  the  eyes  are  viewing  a  3-dimensional  scene.  As  before,  the  disparity  is  the 
difference  between  the  angular  separations  ©1  and  ©2.  For  the  typical  case,  where  ©1  and  ©2  arc 
small. 


Ai  =  Dsin©]  =  D©i 
A2  *  Dsin©2  =  D©2 

Ti  =  ©1  -  ©2  =  ”  (radians)  (5) 
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From  equation  (5),  then,  it  is  clear  that  the  disparity  of  a  point  is  directly  proponional  to  the 
difference  in  its  horizontal  position  in  the  two  images.  Therefore,  as  long  as  the  disparity  does  not 
exceed  the  detectable  range,  large  horizontal  shifts  of  the  target  should  be  easier  to  detect  than 
smaller  shifts. 

Stereoscopes 

When  the  eyes  observe  two  dissimilar  scenes  as  in  Figure  19  muscular  movements  must  be 
to  correctly  align  the  images  on  the  retinae.  If  the  scenes  are  large  and  cannot  be  placed  very 
close  to  one  another,  these  movements  are  very  difficult,  and  correct  alignment  of  the  images  may 
be  impossible  to  obtain.  The  introduction  of  the  first  stereoscope  in  1838  by  Sir  Charles 
Wheatstone  alleviated  this  problem.  Many  stereoscopes  have  been  developed  since  then,  all  having 
the  same  purpose.  The  stereoscope  used  at  the  Valley  Forge  Research  Center  is  based  on  a  design 
proposed  by  Sir  David  Brewster  in  1849  and  later  improved  by  Oliver  Wendell  Holmes.  It  is 
comprised  of  a  pair  of  lenses  and  its  function  is  demonstrated  in  Figure  4. 

In  Figure  20,  the  normal  vision  of  a  single  photograph  is  shown.  Remember  that  when 
both  eyes  see  the  same  image,  as  they  do  here,  there  are  no  disparities  and,  as  a  result,  no  depth  is 
perceived.  In  Figure  21,  we  demonstrate  the  use  of  the  stereoscope  to  view  a  stereo  pair.  The 
image  of  each  photograph,  denoted  by  the  solid  lines,  is  refracted  at  the  lenses  before  entering  the 
eyes.  Notice  that  the  positioning  of  the  eyes  and  the  retinal  placement  of  the  images  are  exactly  the 
same  as  in  Figure  20.  As  a  result,  the  observer  believes  he  is  looking  at  a  single  photograph 
outlined  by  the  dashed  lines.  Now,  however,  since  the  images  originate  from  dissimilar 
photographs,  disparities  are  present  and  the  observer  perceives  the  photograph  in  depth. 

In  working  with  the  stereoscope,  it  has  been  observed  that  the  separation  of  the  target  and 
clutter  may  sometimes  be  more  easily  detected  by  looking  into  the  lenses  at  an  angle,  rather  than 
head  on.  Another  observation  is  that  it  takes  some  people  considerable  amounts  of  time  in  their 
fu^t  attempts  to  fuse  the  images.  But  after  their  eyes  have  "learned"  to  correctly  align  the  images, 
their  ability  improves.  These  findings  suggest  that  an  individual's  ability  to  successfully  detect 
targets  in  clutter  may  improve  with  training. 

POTENTIAL  APPUCATIONS  OF  3-D  RADAR  IMAGING 
Detection  Of  Moving  Targets 

The  experiments  discussed  emulate  this  scenario.  A  single  radio  camera  is  mounted  on  a 
satellite  or  high  flying  aircraft  and  aimed  towards  the  ground.  Successive  pictures  are  taken  and 
observed  stereoscopically.  The  radar  returns  fix)m  the  ground  will  be  identical  in  the  two  images 
and  therefore  will  not  be  observed  in  depth.  However  if  a  target  has  moved  during  the  time  in 
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which  the  pictures  are  taken,  it  will  appear  in  different  locations  in  the  two  images  and  will  be 
observed  separate  from  the  clutter. 

When  the  movement  of  the  target  can  be  displayed  as  a  purely  horizontal  shift,  as  in  the 
experiment,  the  size  of  the  shift  is  proportional  to  the  speed  of  the  target.  Since  the  perceived 
distance  between  the  target  and  clutter  is  proportional  to  this  shift,  the  speed  of  the  target  can  be 
estimated  from  the  seperation  of  the  target  and  clutter  in  the3-D  image.  In  this  manner,  targets  with 
«;mall  radar  cToss-sections  can  be  differentiated  by  the  speeds  at  which  they  are  moving. 

Detection  of  Static  Targets 

Two  radio  cameras  focused  on  the  same  location  simulate  eyes  operating  at  microwave 
rather  than  optical  frequencies.  In  applications  such  as  locating  targets  below  the  surfaces  of  land 
and  water,  the  target  can  be  located  in  both  position  and  depth.  In  another  scenario,  two  cameras 
mounted  to  satellites  could  determine  the  locations  of  flying  aircraft  since  the  aircraft  targets  would 
be  observed  above  the  earth  background.. 


CONCLUSIONS 

It  is  clear  that  there  is  great  potential  for  three-dimensional  radar  systems  in  the  area  of 
detection.  By  increasing  the  dimensionality  of  the  detection  system,  it  has  become  possible  to 
identify  targets  that  would  not  be  detected  otherwise.  During  this  period,  we  have  isolated  the 
factors  necessary  to  obtain  this  high  level  of  detection,  namely  high-resolution  and  large  disparity. 


W)  (CORRECTED) 
(M)(VmCQPa58) 
2-0  (VFRC  OPR.  58) 


Figure  17 , 2-D  detection  versus  0-D  as  a  function  of  TCR. 
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Figure  l?  Stereo  Vision  of  a 
Three-Dimensional  Scene. 


Figure  19  Reproduction  of  a 
iWe-Dimensional  Scene. 


Figure  20 
Normal  Vision  of 
a  Photograph. 


Vision  of  a  Stereo  Pair  With  Lenses. 
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5. 


AUTOMATIC  TARGET  RECOGI^TTION  USING  NEURAL  NETWORKS  FOR 
MICROWAVE  IMAGING  SYSTEMS;  PRELIMINARY  REPORT 


OBJECTIVE 

Two-dimensional,  high  resolution  microwave  imagery  has  proven  to  be  a  powerful  tool  in  the 
recognition  of  radar  targets  [1].  The  objective  of  the  proposed  effort  is  to  extend  this  recognitive 
capability  to  automatic  target  recognition  by  using  artificial  neural  networks.  The  neural  networks 
developed  will  be  able  to  correctly  recognize  and  classify  fragmented  (incomplete),  distorted,  or  noise- 
corrupted  imagery  obtained  with  microwave  imaging  systems. 

INTKODUCnON 

Microwave  imaging  systems,  which  combine  radar  with  adaptive  signal  processing  techniques, 
have  achieved  two-dimensional  high  resolution  images  of  aircraft.  The  processing  is  based  on 
advanced  signal -imaging  algorithms  [  i  J.  Several  radio  camera  systems  currently  in  operation  at  the 
Valley  Forge  Research  Center  of  the  Moore  School  demonstrate  the  capability  of  microwave  imaging 
systems  to  produce  images  suitable  for  target  recognition.  Fig.  22  illustrates  a  high  quality  two- 
dimensional  image  of  an  L- 10 11  approaching  the  Philadelphia  International  Airport.  A  plan-view 
drawing  of  the  airplane  is  shown  in  Fig.  22  to  facilitate  comparison. 


Figure  22  Microwave  image  of  L-101 1  obtained  with  microwave  imaging  system:  (A)  Microwave 
image.  (B)  Plan-view.  The  radio  camera  operated  at  X-band  (3  cm  wavelength)  with 
150  MHz  bandwidth.  The  aircraft  distance  was  approximately  2.5  KM. 

The  high  quality  exemplified  by  the  microwave  image  of  Fig.  22  however,  is  sometimes 
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lacking  because  of  factors  such  as  random  glint,  scintillation,  background  clutter,  inadequate  dynamic 
range,  system  noise,  and  compensation  error.  The  images  may  be  fragmented,  distorted,  or  severely 
corrupted  by  noise,  to  the  extent  that  a  human  observer  is  unable  to  recognize  the  target  Two 
examples  are  shown  in  Fig.  23 . 
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Figure  23 ,  Two  examples  of  fragmented  microwave  images :  an  L-  101 1  and  a  BOEING  727. 


In  order  to  improve  the  2-D  images,  diversity  combining  techniques  may  be  introduced  in  the 
imaging  procedure.  Several  incomplete  “raw”  images  derived  from  different  angles,  frequencies,  or 
polarizations  can  form  a  superior  image  by  overlapping  or  averaging.  Fig.24shows  the  diversity 
process  exercised  on  several  images  of  a  Boeing  727.  The  upper  row  shows  five,  independent  “raw” 
images,  each  highly  fragmented.  The  lower  row  shows  that  2nd  and  5th  order  diversity  lead  to  less 
fragmented  images.  Unfortunately,  this  technique  is  not  practical  for  those  applications  in  which  the 
time  required  for  data  acquisition  and  processing  is  too  long.  This  problem  becomes  even  more 
significant  when  multiple  frequencies  and  polarizations  arc  introduced.  Thus  the  ability  of  neural  nets 
to  recognize  highly  fragmented  targets  is  very  important 

In  theory,  obtaining  a  perfect  2-D  image  of  a  target  requires  an  infinite  antenna  array  to  coUea 
all  the  information  over  the  entire  range  of  frequency  and  polarization.  Using  an  infinite  array  means 
that  a  5-function  radiation  pattern  scans  across  the  target-  Measuring  over  the  entire  frequency 
range  means  that  the  target  is  illuminated  by  an  impulse  of  infinitesimal  extent  to  obtain  extremely  high 
range  resolution.  Measuring  with  all  polarizations  produces  still  more  information  about  the  target  In 
practical  applications,  however,  the  amount  of  information  that  can  be  gathered  is  always  limited.  A 
2-D  image  is  formed  only  from  partially  available  information  in  a  noisy  environment 


In  addition,  the  2-D  images  obtained  with  microwave  imaging  systems  are  only  projections  in 
range  and  angle  of  the  original  3-D  objects.  Different  projections  will  be  produced  at  different 
elevations.  In  other  words,  the  image  of  an  object  formed  directly  with  a  microwave  imaging  system 
is  a  deformed  version  of  the  2-D  plan  view  of  the  object 


Figure  2  4  Diversity  process  exercised  on  several  images  of  BOEING  727,  showing  reduced 
fragmentation 


It  is  iiiqxjrtant  for  us  to  resolve  the  problem  of  automadc  recognition  with  limited,  incomplete, 
distorted,  or  noise-corrupted  images.  This  is  especially  true  when  the  images  are  not  useful  for 
reliable  human  identification. 

Human  nervous  systems  are  well  known  to  be  good  at  recognidon  fiom  partial,  distorted  and 
obscured  informadon:  They  are  able  to  recall  an  item  from  memory  easily  with  only  partial  informadon 
on  that  item  (  associanve  or  contant-addressabie  memory. )  Artificial  neural  networks  mimic  human 
nervous  systems  based  on  neurophysiological  and  neuroanaiomical  studies:  These  networks  have 
many  of  the  attractive  features  of  the  human  nervous  system  such  as  rapid  processing,  robustness, 
fault  tolerance,  capabilities  to  adapt  and  continue  learning,  and  they  can  perform  many  intelligent 
information-processing  funaons 

So  far,  many  neural  network  models  have  been  developed  and  successfully  applied  to  image 
recognition.  Our  preliminary  investigations  also  demonstrate  that  the  recognition  of  poor 
microwave  imagery  based  on  neural  network  models  is  feasible  and  efficient 
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The  necx:ognitron  was  fashioned  by  analogy  to  visual  systems,  and  therefore  it  could  be 
applied  directly  to  the  recognition  of  low  quality  2-D  microwave  images.  In  the  case  of  airplanes,  the 
prominent  features  include  the  wingspan,  the  length  and  width  of  the  fuselage,  the  locations  and 
numbers  of  engines,  the  orientation  of  the  wings  with  respect  to  the  fuselage,  and  the  shapes  and 
positions  of  the  wings  and  tail  structure. 

One  of  the  proposed  efforts  is  to  make  a  detailed  comparison  of  the  above  three  types  of  neural 
network  models  and  their  relative  applicability  to  the  recognition  of  microwave  images.  The  most 
suitable  naodel  of  the  three  types  will  then  be  selected,  and  improved  if  necessary. 

HARDWARE  IMPLEMENTATION  OF  THE  NEURAL  NETWORK  MODEL 

There  are  three  ways  to  make  practical  implementations  of  neural  networks:  (1)  Develop 
simulation  software  designed  on  a  general-purpose  digital  computer  to  perform  the  actions  of  a  true 
neural  netwoik;  (2)  build  a  special-puipose  digital  neural  network  emulator  that  connects  to  a  general- 
purpose  host  system;  or  (3)  construct  a  true  electronic  or  optical  neural  network.  We  plan  to  build  a 
special-purpose  digital  hardware  emulator  to  realize  our  neural  network  model  because  the  emulator  is 
effective,  low  in  cost  and  easier  to  implement. 

For  this  purpose,  some  advanced  digital  signal  processing  microprocessors,  such  as  the 
ADSP21000  (Analog  Devices,  Inc.),  DSP32C  (AT&T  BeU  Lab.),  DPS96002  (Motorola,  Inc.),  and 
TMS320C30  (Texas  Instruments,  Inc.),  will  be  considered.  They  have  many  of  the  characteristics 
needed  by  a  neural  network  emulator.  For  example,  the  TMS 320030  has  a  very  high  speed  ( 60  ns 
cycle-time ).  Its  very  efficient  pipeline  function  can  execute  each  instruction  in  a  single  cycle.  Its 
floating-point  arithmetic  capability  permits  the  handling  of  numbers  of  very  high  dynamic  range.  Its 
parallel  (  single  cycle  )  multiplier  instructions  and  block-repeat  capability  will  be  very  useful  in  the 
implementation  of  the  neural  network  emulator.  In  addition,  it  has  one  4K  x  32  bits  ROM  block 
and  two  IK  x  32  bits  RAM  blocks.  The  separate  program  buses,  data  buses,  and  DMA  buses  allow 
for  parallel  program  fetches,  data  reads  and  writes,  and  DMA  operations.  Because  of  the  on-chip 
DMA  controller,  the  TMS320C30  can  interface  to  slow  external  memories  and  peripherals  without 
reducing  throughput  to  the  CPU.  In  multilayer  neural  networks,  each  layer  can  use  a  chip  to 
implement  its  actions.  Using  this  type  of  microprocessor,  we  will  develop  an  advanced  neural  network 
emulator  to  implement  the  automatic  recognition  of  images  from  microwave  imaging  systems. 

RESULT  OF  PRELIMINARY  INVESTIGATIONS 

We  used  simulated  data  sets  and  an  experimental  data  base  of  high  resolution  microwave 
imaaes  of  commercial  aircraft  (Boeing  727,  Boeing  747,  L-101 1  and  DC-8)  flying  into  Philadelphia 
International  Airpon.  The  radar  operated  at  X-band  (k=3  cm)  with  150  MHz  bandwidth.  The 


resolution  cell  was  approximately  1  m^.  The  target  distances  from  the  Valley  Forge  radio  camera  are 
approximately  3  km- 

We  tested  the  Hamming  net  and  the  multilayer  pcrceptron  in  our  investigations.  The  airplane 
images  were  composed  of  59  x  70  pixels.  When  the  Hamming  net  was  used,  the  images  were  first 
transformed  to  binary  format .  Their  exemplars  were  plan-view  drawings  obtained  from  an  airplane 
manual.  The  inputs  were  fragmented  microwave  images  to  which  noise  has  been  added.  Fig  25 
shows  typical  results.  The  first  line  in  Fig.25  shows  the  fragmented  images  without  noise  and  the 
second  line  shows  them  after  noise  has  been  added.  Noise  was  added  by  randomly  changing  white 
pixels  to  black  and  black  pixels  to  white  with  a  0.3  probability.  This  is  equivalent  to  a  SNR  of  3.7 
dB.  The  third  line  shows  the  exemplars  chosen  by  the  Hamming  network  that  best  match  the  noisy 
and  fragmented  aircraft  images.  In  each  case  the  selection  was  the  correct  choice  out  of  four. 
Experiments  consistently  demonstrated  that  the  Hamming  net  can  work  well  if  the  input  imagery  is 
identical  in  size,  position  and  orientation  to  its  exemplar.  If  the  input  images  have  minor  changes  in 
size,  position  and  orientation,  the  net  still  can  recognize  them  correctly.  When  these  changes  are 
large,  however,  mistakes  result.  The  effects  of  variations  in  the  scale,  position  and  orientation  is  an 
important  matter  to  be  studied. 

Unlike  the  Hamming  net,  the  Perceptron  requires  training.  When  the  three-layer  perceptron 
with  a  rotation-independent  net  was  used,  the  training  set  used  in  the  experiment  had  20  images  of 
these  commercial  aircraft  These  included  two  real  images  produced  in  the  Valley  Forge  Research 
Center  and  three  simulated  images  of  each  type.  The  simulated  training  images  were  pieces  of  the 
exemplars,  with  up  to  50%  random  fragmentation.  The  neural  network  was  trained  using  the 
backpropagation  learning  algorithm  [  14  ].  Internal  thresholds  were  set  to  equal  0.  During  the 
training  process  the  weight  values  were  changed  to  improve  performance.  The  trained  network  was 
then  presented  with  eight  highly  fragmented  (but  noise  free)  test  images  excluded  from  the  training  data 
sets.  Fig.  26  shows  the  eight  test  images.  The  first  line  in  Fig.  2  6  shows  four  real  images  and  the 
second  line  shows  four,  simulated  images.  The  test  results  showed  that  this  network  has  good  ability 
to  generalize,  since  each  of  the  test  images  was  successfully  detected  and  recognized.  The  three-layer 
perceptron  with  a  rotation-independent  net  can  learn  to  recognize  imagery  in  the  training  data  set,  and 
the  trained  network  could  be  used  to  correctly  recognize  imagery  in  the  test  data  set 
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(B)  SIMULATED  TEST  IMAGES 
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